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Abstract 



The measurement of several event shape variables is presented in minimum bias pp collisions 
at a center-of-mass energy of 7 TeV using the ATLAS detector at the LHC. The observables studied 
are the transverse thrust, thrust minor and transverse sphericity each defined using the momenta 
perpendicular to the beam direction of the final state charged particles. Events with at least six 
charged particles are selected. In addition to the differential distributions, the evolution of each event 
shape variable as a function of the leading charged particle transverse momentum, charged particle 
multiplicity and summed transverse momentum is presented. Predictions from several Monte Carlo 
models have been compared with data and they show significant deviations from data. 
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I. INTRODUCTION 

Event shape variables describe the properties of the 
energy flow and the structure of hadronic events. In this 
analysis, three event shape observables [1, 2] are mea- 
sured: the transverse thrust, the thrust minor and the 
transverse sphericity, each built from the charged particle 
momenta using tracking information from proton-proton 
collisions at i/s = 7 TeV collected with the ATLAS detec- 
tor [3]. Event shape observables are among the simplest 
experimentally measured quantities and depending on 
the events being considered, may have sensitivity to both 
perturbative and non-perturbative aspects of Quantum 
Chromodynamics (QCD). 

Event shapes have been investigated extensively in 
e + e~ and ep deep-inelastic scattering experiments to 
study the energy flow in the hadronic final state, test 
the predictions of perturbative QCD and to extract a pre- 
cise value for the strong coupling constant a s [4-11]. In 
pp collisions at the Tevatron, the dependence of the event 
shape observables on the transverse energy of the leading 
jet and on contributions from the underlying event has 
been studied [12]. At the Large Hadron Collider (LHC), 
event shape observables have recently been studied in 
inclusive interactions [13] and multijet events [14, 15]. 

The study of event shape observables in inclusive 
minimum-bias events plays an important role in under- 
standing the nature of soft-QCD processes at the LHC 
center-of-mass energies [16], where "soft" refers to interac- 
tions with low momentum transfer between the scattering 
particles. Soft interactions cannot be reliably calculated 
from theory and are thus generally described by phe- 
nomenological models, usually implemented in Monte 
Carlo (MC) event generators. These models contain 
many parameters whose values are a priori unknown 
and thus need to be constrained by measurements. Inclu- 
sive and semi-inclusive observables sensitive to soft-QCD 
processes have been measured at the LHC by ATLAS [17- 
19], CMS [20, 21] and ALICE [22, 23]. The measurements 
presented in this paper will further constrain the event 
generator models, which encapsulate our understanding 
of these soft processes. 



In this analysis, the event shape observables are con- 
structed from six or more primary charged particles in 
the pseudorapidity range \r)\ < 2.5 and with transverse 
momentum px > 0.5 GeV [24]. Primary charged par- 
ticles are defined as those with a mean proper lifetime 
r > 30 ps, produced either directly in the pp interaction 
or from the subsequent decay of particles with a shorter 
lifetime. The particle level refers to particles as they 
emerge from the proton-proton interaction. The detector 
level corresponds to tracks as measured by the detector(s) 
after interaction with the material, and including the 
detector response. The results are corrected for detector 
effects, using simulation, to obtain distributions of the 
event shape variables defined at particle level which can 
be directly compared to MC models. 

This paper is organized as follows: Section II defines 
the event shape variables; the detector is described in 
Section III; Section IV discusses the MC models used in 
this analysis; Section V and Section VI respectively de- 
scribe the event selections and background contributions. 
The correction of the data back to particle level, and 
estimation of the systematic uncertainties are described 
in Section VII and Section VIII; the results are discussed 
in Section IX and finally the conclusions are presented in 
Section X. 



II. EVENT SHAPE OBSERVABLES 

In particle collisions, event shape observables describe 
the geometric properties of the energy flow in the final 
state. A single event shape variable can distinguish in a 
continuous way between configurations in which all the 
particles are flowing (forward and backward) along a single 
axis and configurations where the energy is distributed 
more uniformly over the 47r solid angle. If defined as the 
ratio of measured quantities, the corresponding systematic 
uncertainties may be small. 

In hadron collisions, where the center-of-mass frame of 
the hard interaction is usually boosted along the beam 
axis, event shape observables are often defined in terms 
of the transverse momenta, which are Lorentz-invariant 
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under such boosts. Different formulations of event shape 
observables are possible; the most intuitive is to calculate 
the event shape from all particles in an event. These are 
denoted as directly global event shapes [1, 2]. In most 
hadronic collisions, it is not usually possible to detect 
all particles in an event due to the finite detector accep- 
tance limited at small angles by the presence of the beam 
pipe. Event shapes which include only particles from 
a restricted phase space in pseudorapidity n, are called 
central event shapes: in this analysis charged particles 
within 1 77 1 < 2.5 are used. These central event shape val- 
ues are nevertheless sensitive to non-perturbative effects 
at low momentum transfer and can be experimentally 
useful observables. 

The thrust is one of the most widely used event shape 
variables. At hadron colliders the transverse thrust for a 
given event is defined as: 



T± = max 



i 



(1) 



where the sum is performed over the transverse momenta 
pr,i of all particles in the event. The thrust axis fix is 
the unit vector n that maximizes the ratio in Eq. 1. The 
transverse thrust ranges from T± = 1 for a perfectly bal- 
anced, pencil-like, dijet topology to T± = (| cos-0|) = 2/tt 
for a circularly symmetric distribution of particles in the 
transverse plane, where ip is the azimuthal angle between 
the thrust axis and each respective particle. Since many 
event shape variables vanish in a balanced dijet topol- 
ogy, it is convenient to define tj_ = 1 — Tj_ , which shares 
this property. Hereafter, any discussion of the observable 
called transverse thrust will refer to the quantity rj_. 

The thrust axis nx and the beam axis z define the 
event plane. The transverse thrust minor measures the 
out-of-event-plane energy flow: 

£ lFr,i-n m | 

T M = ' - , n m = n T x z . (2) 

2jPT.il 



values of this tensor, A2 and A3 
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The sphericity has values between and 1, where a bal- 
anced dijet event corresponds to S — and an isotropic 
event to S = 1. Sphericity is essentially a measure of 
the summed p\ with respect to the event axis [25, 26], 
where the event axis is defined as the line passing through 
the interaction point and oriented along the eigenvector 
associated with the largest eigenvalue Ai. Similar to trans- 
verse thrust, the transverse sphericity, S±, is defined in 
terms of the transverse components only: 



S xy = £ 



1 
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where A^ < A* y are the two eigenvalues of S xy . 
The following distributions are measured: 



• Normalized distributions: 



dN c , 



1 dN ev 

W„ dr<i h ' JV CV dT^ 1 



• Average values: (rj ), (TjJ) and (Sjf) as functions 
of N ch and J^Pt'i 

where iV cv is the number of events with six or more 
charged particles within the selected kinematic range; 
iVch is the number of charged particles in an event; J^Pt 
is the scalar sum of the transverse momenta of the charged 
particles in the event. The event shape observables rf 1 , 
T^J 1 and Sf 1 are defined as above, with the superscript 
indicating that they are constructed from charged parti- 
cles. The three normalized differential distributions are 
studied separately for p!p ad above 0.5 GeV, 2.5 GeV and 
5 GeV, where pi£ ad is the transverse momentum of the 
highest pt (leading) charged particle. 



III. THE ATLAS DETECTOR 



The transverse thrust minor is for a pencil-like event in 
4> and 2/-7T for an isotropic event. 

Another widely used event shape variable is the spheric- 
ity, S, which describes the event energy flow based on the 
momentum tensor, 



S 



af3 _ 



(3) 



where the Greek indices represent the x, y, and z compo- 
nents of the momentum of the particle i. The sphericity 
of the event is defined in terms of the two smallest eigen- 



The ATLAS detector [3] covers almost the full solid 
angle around the collision point with layers of tracking 
detectors, calorimeters and muon chambers. The compo- 
nents that are relevant for this analysis are the tracking 
detectors. The inner tracking detector has full coverage 
in azimuthal angle cj> and covers the pseudorapidity range 
\n\ < 2.5. It consists of a silicon pixel detector (pixel), 
a semiconductor tracker (SCT) and a straw-tube transi- 
tion radiation tracker (TRT). These detectors are located 
at a radial distance from the beam line of 50.5-150 mm, 
299-560 mm and 563-1066 mm, respectively, and are im- 
mersed in a 2 T axial magnetic field. The inner detec- 
tor barrel (end-caps) consist of 3 (2 x 3) pixel layers, 4 
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(2 x 9) layers of double-sided silicon strip modules, and 
73 (2 x 160) layers of TRT straw-tubes. These detectors 
have position resolutions typically of 10 urn, 17 um and 
130 pm for the r-<j> coordinate. The pixel and SCT de- 
tectors provide measurements of the r-z coordinate with 
typical resolutions of 115 pm and 580 pm, respectively. 
The TRT acceptance is \r)\ < 2.0. A track traversing the 
barrel typically has 11 silicon hits (3 pixel clusters and 8 
strip clusters) and more than 30 straw-tube hits. 

The measurements presented here rely on the minimum- 
bias trigger scintillator (MBTS) system. The MBTS 
detectors are mounted at each end of the tracking detector 
at z — ±3.56 m and are segmented into eight sectors 
in azimuth and two concentric rings in pseudorapidity 
(2.09 < M < 2.82 and 2.82 < \r\\ < 3.84). The MBTS 
trigger was configured to require at least one hit above 
threshold from either side of the detector in coincidence 
with a fast beam pickup device ensuring that the event is 
compatible with a bunch crossing. 



IV. MONTE CARLO MODELS 

Monte Carlo event samples are used to compute detec- 
tor acceptance and reconstruction efficiency, determine 
background contributions, correct the measurements for 
detector effects, and to calculate systematic uncertainties. 
Finally, different phenomenological models implemented 
in the MC generators are compared to the data corrected 
to the particle level. 

The PYTHIA6 [27], pythia8 [28] and Herwig++ [29, 30] 
event generators were used to produce the simulated event 
samples for the analysis. These generators implement 
leading-logarithmic parton shower models matched to 
leading-order matrix element calculations with different 
hadronization models and ordering for the parton shower. 
The PYTHIA6 and PYTHIA8 generators use a hadronization 
model based upon fragmentation of color strings and 
a px-ordered or virtuality ordered shower, whereas the 
Herwig+- 1- generator implements a cluster hadronization 
scheme with parton showering ordered by emission angle. 
Different settings of model parameters optimized, or tuned 
to reproduce the existing experimental data have been 
used for the MC generators. 

The following Monte Carlo generators were used: 
pythia6 with tunes AMBT1 [31], AMBT2B [32], DW [33], 
Zl [34]; pythia 8 tune A2 [35]; and Herwig++ tune UE7- 
2 [36]. For PYTHIA 6, version 6.425 was used in all cases 
except for tune DW, for which version 6.421 was used. 

The AMBT1 (ATLAS Minimum Bias Tune 1) is the 
first LHC data tune of pythia 6 from ATLAS, and uses 
the diffraction suppressed part of the early minimum-bias 
measurements [17]. This employs the MRST LO** [37] 
parton distribution functions (PDFs) and the pythia 6 
Px-ordered parton shower. The pythia 6 DW tune uses 
a virtuality ordered parton shower and a multi-parton 
interaction (MPI) model not interleaved with the initial 
state radiation (ISR). This tune was constructed to de- 



scribe CDF Run II underlying event data, and uses the 
CTEQ5L1 PDF set. The AMBT2B tune is an improved 
version of AMBT1 tune, which optimizes parameters con- 
trolling the ISR cutoff and evolution by including ATLAS 
track-jet [38], jet shape [39] and dijet decorrelation [40] 
data. The tune with CTEQ6L1 PDF is used here, and 
provides one of the best available descriptions of existing 
minimum-bias data at y/s = 7 TeV. The Zl tune devel- 
oped by the CMS collaboration with the CTEQ5L PDF 
set is based on AMBT1, but uses CMS charged-particle 
jet underlying event data to obtain a better description 
of the underlying event. The pythia 8 generator uses the 
MPI model interleaved with both ISR and final state radi- 
ation (FSR). The ATLAS minimum-bias tune of pythia 8, 
A2 with the MSTW2008LO PDF, which provides a good 
description of minimum-bias data at y^s = 7 TeV, is used. 
For Herwig++, version 2.5.1 is used with a 7 TeV underly- 
ing event tune, UE7-2, which employs color reconnection, 
and uses the MRST LO** PDF. 

The reference tune for this analysis is chosen to be 
pythia 6 AMBT1. Samples generated with this tune 
were passed through the ATLAS detector and trigger sim- 
ulations [41] and then reconstructed and analyzed using 
the same procedure and software that are used for the 
data. Reconstructed MC events are then used to correct 
the data for detector effects. A sample generated with an 
older version of Herwig++, 2.5.0 with no additional tun- 
ing, was also passed through the full detector simulation 
and the analysis chain for systematic studies of unfolding 
corrections. 



V. EVENT AND TRACK SELECTION 

The collision data used for the analysis presented here 
were collected in April 2010 with a minimal prescale factor 
for the minimum-bias trigger. The only further require- 
ment for selecting the data sample is that the MBTS 
trigger and all inner detector subsystems were at nominal 
operating conditions. In each event the reconstructed 
vertices are ordered by the p\ over the tracks assigned 
to each vertex, and the vertex with the highest ^ p\ is 
taken as the primary interaction vertex of the event. To 
reduce the contribution from beam-related backgrounds 
and decays of long-lived particles, and to minimize the 
systematic uncertainties, events are rejected if they con- 
tain any other vertex reconstructed with four or more 
tracks. 

If there is only one vertex in the event, or if any addi- 
tional vertex in the event has three or fewer tracks, all 
tracks from the event that pass the track selection de- 
scribed next are retained. After this selection, the fraction 
of events with more than one proton-proton interaction 
in the same bunch crossing (referred to as pile-up) is 
found to be approximately 0.1% and this residual contri- 
bution is therefore neglected. The average number of pp 
interactions per bunch crossing during this data taking 
period was less than 0.15, indicating a negligible pile- 
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up contribution. The MC samples used have no pile-up 
contribution. 

Events are required to contain at least six tracks that 
fulfill the following criteria: 

• p T > 0.5 GeV; 

• a minimum of one pixel and six SCT hits; 

• a hit in the innermost pixel layer, if the correspond- 
ing pixel module was active; 

• transverse and longitudinal impact parameters with 
respect to the primary vertex, |d | < 1.5 mm and 
\zq\ sin 6* < 1.5 mm; 

• a track- fit probability \ 2 > 0.01 in order to remove 
mis-measured tracks for tracks with px > 10 GeV. 

Tracks with p^ > 0.5 GeV are less prone than lower- 
Pt tracks to inefficiencies and systematic uncertainties 
resulting from interactions with the material inside the 
tracking volume. 

After event selection, the analysis is based on approx- 
imately 17 million events containing approximately 300 
million tracks. For the PYTHIA 6 generator and for the 
pythia 8 generator, which has a harder diffractive model 
than the former, the contribution to the event shape 
observables from diffractive events is negligible when re- 
quiring six or more tracks in the event. 



VI. BACKGROUND CONTRIBUTIONS 

A. Backgrounds 

Backgrounds comprise beam-induced events, due to 
beam-gas and beam-material interactions, as well as non- 
beam background, from cosmic-ray interactions and de- 
tector noise. The contribution of these background events 
remaining after the event selection is estimated using the 
number of pixel hits not associated with reconstructed 
tracks. This multiplicity includes unassigned hits from 
1ow-pt looping tracks, but is dominated at higher multi- 
plicities by hits from charged particles produced in inelas- 
tic interactions of protons with the residual gas inside the 
beam pipe. The vertex requirement removes most of the 
beam background events and the residual contribution is 
below 0.1%. As the level of background is very low, no 
explicit background subtraction was performed. 



B. Secondary track fraction 

The primary charged particle multiplicities are mea- 
sured from selected tracks after correcting for the fractions 
of secondary and poorly-reconstructed tracks in the sam- 
ple. The potential background from fake tracks is found 
from MC studies to be less than 0.01% [17]. 



Non-primary tracks arise predominantly from hadronic 
interactions, photon conversions to positron-electron pairs 
in the detector material and decays of long-lived particles. 
For > 0.5 GeV the contribution from photon conver- 
sions is small. The systematic uncertainty from secondary 
decays is included in the uncertainties associated with the 
tracking performance. 

VII. CORRECTION TO PARTICLE LEVEL 

To facilitate comparison with theoretical predictions 
and other measurements, the event shape distributions 
are presented at particle level for charged particles, after 
correction for trigger and event selection efficiencies, as 
well as detector resolution effects. A two-step correction 
procedure is used: first, corrections for event efficiency are 
applied, followed by an additional bin-by-bin correction to 
account for tracking inefficiencies, possible bin migrations 
and any remaining detector effects. 

A. Event-level correction 

Trigger and vertexing efficiencies are taken from a pre- 
vious analysis using the same data sample [17]. The 
efficiency of the MBTS trigger is determined from data 
using a control trigger and found to be fully efficient for 
the analysis requirement of at least six tracks. The ver- 
tex reconstruction efficiency is also measured in data by 
taking the ratio of the number of triggered events with 
a reconstructed vertex to the total number of triggered 
events. This ratio is also found to be very close to unity. 
The total correction applied to account for events lost 
due to the trigger and vertex requirements is less than 
1% and it varies very weakly with the number of tracks 
associated with the primary vertex. 

B. Bin-by-bin correction 

The event shape observables presented here are sensitive 
to changes in the configuration of the selected tracks. 
Applying average track efficiencies to individual tracks 
on a track-by-track basis and reweighting tracks distorts 
the event shape distribution. A more robust approach is 
to apply bin-by-bin corrections to find the event shape 
distribution at particle level. Such a bin-by-bin correction 
is applied to all distributions after applying the event-level 
efficiency corrections described above. 

The correction factors Cbi n are evaluated separately in 
each bin for each event shape observable, 

j/Gen 

ft bin rj\ 

bln T/ Rcco, off coir ' w 

where V^° n and V^ co ' off corr represent the generator level 
MC value of the bin content and the reconstructed MC 
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value after applying the event-level efficiency corrections 
for each bin, respectively. The corrected value of the 
bin content for an observable is found by multiplying the 
measured bin content by the corresponding correction 
factor. The bin sizes are chosen to be consistent with the 
resolution of the correction procedure. 

The correction factors are calculated using the two 
different models implemented in pythia 6 AMBT1 and 
Herwig++. This correction accounts for bin- by-bin mi- 
grations and tracking inefficiencies. For each distribution, 
the unfolding factor is typically within ±f0% of unity 
for most of the range and very close to unity for the 
average values. The difference from unity becomes more 
pronounced at the statistically limited edges of the distri- 
butions. The correction factors for the inclusive distribu- 
tions of the three event shape observables are shown in the 
bottom panels of Fig. f for the two MC event generators 
mentioned above. Though the two MC generators have 
different distributions, the bin- by-bin correction factors 
are similar. 



VIII. SYSTEMATIC UNCERTAINTIES 

Systematic uncertainties on the measured distributions 
are assessed with the following sources of uncertainty 
included: 

Tracking: The largest of the systematic uncertainties for 
the tracking inefficiency [17] are found to be due to 
the material description in the inner detector. This 
is determined to affect the efficiency by a relative 
difference of 2% in the barrel region, rising to ~ 7% 
for 2.3 < \rj\ < 2.5, for tracks with p T > 0.5 GeV. 
The contribution of the propagated uncertainty to 
the distributions of the event shape variables is 
found to be less than 1 % . 

Bin-by-bin correction model dependence: The re- 
maining contributions to the overall systematic un- 
certainty result from the specific correction method 
used in this analysis. The bin-by-bin corrections in 
general depend on the number of charged particles 
and their pt distributions, so there is some depen- 
dence on the event generators. In order to estimate 
this uncertainty, it is necessary to compare different 
plausible event generators, which deviate signifi- 
cantly from each other, but still give predictions 
close to the data. The corrected results using the 
two very different pythia 6 AMBT1 and Herwig++ 
models are compared. As these two generators use 
very different soft-QCD models the difference is as- 
signed as a systematic uncertainty. The generated 
and reconstructed distributions are shown in Fig. f 
for the two MC event generators and compared with 
the detector-level data. Since this uncertainty is in- 
dependent of any efficiency systematic uncertainties, 
it is added in quadrature with the efficiency system- 
atic uncertainty and the statistical uncertainty. 
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FIG. 1. The generated and reconstructed MC distributions for 
transverse thrust, thrust minor and sphericity are shown in 
the top half of each plot. The correction factors are shown at 
the bottom half for pythia 6 AMBT1 and Herwig+4- default 
tune. The data are shown with only the efficiency corrections 
and statistical uncertainties. 
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Statistical uncertainty of bin-by-bin correction: 

In addition to the model-dependent uncertainty in 
the bin-by-bin correction, there is also a statistical 
uncertainty due to the finite size of the MC sample. 
The statistical fluctuations of the pythia6 AMBT1 
correction factor are found to be negligible for most 
of the distributions, increasing to a few percent in 
the tails of the distributions. This is also added to 
the overall systematic uncertainty estimate. 

The systematic uncertainty due to the small number of 
residual multiple proton-proton interactions is estimated 
to be negligible. 

Table I lists representative values for the various contri- 
butions to the systematic uncertainties for all the event 
shape observables away from the edges of the distribu- 
tions. 



TABLE I. Summary of systematic uncertainties in %. 



Trigger and vertex efficiency 


< 


0.1 


Track reconstruction 


0.1 - 


0.5 


Correction model difference 


1 


- 5 


PYTHIA correction stat. uncertainty 


0.1 


- 2 


Total systematic uncertainty 


1 


- 5 



IX. RESULTS AND DISCUSSION 

The distributions of the transverse thrust, thrust minor 
and transverse sphericity are presented in Figs, 2 and 3, for 
three different values of the minimum p!^ ad . As in previous 
measurements [18], the p T of the leading charged particle 
is used to probe the emergence of a jet-like structure. The 
behavior of the average values of the shape variables as 
functions of the charged particle multiplicity, A c h, and 
transverse momentum scalar sum, ^pt, is presented in 
Fig. 4. Predictions from pythia 6 AMBT2B, pythia 
6 DW, pythia 6 Zl, pythia 8 A2 and Herwig++ UE7- 
2 models are also shown. AMBT2 is chosen instead of 
AMBT1 used to correct the data back to the particle level 
because it shows a slight improvement in reproducing the 
distributions of charged particle transverse momentum 
and multiplicity. 

The inclusive distributions shown in Figs. 2 and 3 indi- 
cate a prevalence of high sphericity events. The shape of 
the transverse thrust and thrust minor distributions does 
not change appreciably upon increasing the lower limit of 
p!j3 ad to 2.5 GeV, while a slight shift toward less spherical 
events and a broadening of the distributions is observed 
for events with pig ad > 5 GeV. The mean of the transverse 
thrust and thrust minor respectively changes from 0.23 
and 0.51 for inclusive distributions to 0.22 and 0.48 for 
distributions with pJp d > 5 GeV. The transverse thrust 



and thrust minor have little sensitivity to the leading 
charged particle transverse momentum threshold in the 
explored range. In contrast, the distribution of transverse 
sphericity shows a clear transition from spherical toward 
dijet-like events in the same pl£ ad range. 

The mean values of event shape observables as functions 
of Af c h and YIpt are shown in Fig. 4. They are seen to 
increase to values largely consistent with the position of 
the maxima of the corresponding differential distributions 
with pig ad > 0.5 GeV. Before reaching the plateau, for 
low values of A c h and J^pt, the mean values of the event 
shape variables correspond to less spherical events. This 
can also be seen in the differential distributions as the 
lower limit on the leading charged particle transverse 
momentum is reduced. At very high values of A c h and 
J^Pt, the mean values indicate more spherical events. 

In general, the MC models predict fewer high sphericity 
events than the data. However, the Zl tune does a rea- 
sonable job of reproducing all of the distributions. The 
AMBT2B prediction shows better agreement for the in- 
clusive distributions in Figs. 2(a), 2(b) and 3(a), than the 
PYTHIA 8 A2 tune and Herwig++ UE7-2 tune predictions. 
The pythia 6 DW tune predictions are consistently fur- 
thest from the data, as observed for the charged particle 
multiplicity and px distributions [17] also. It is interest- 
ing to note that AMBT2B does not describe the event 
shapes very well, even though it has been constructed to 
model the data using inclusive distributions of charged 
particle multiplicity and transverse momentum presented 
in Ref. [17]. This suggests that event shape observables 
may be useful for tuning MC in the future. 

For events with a higher pJp d , the Zl tune provides 
the best description of the data, followed by Herwig++. 
In Figs. 2(c) and 2(d), with p!^ ad > 2.5 GeV, for trans- 
verse thrust and thrust minor, the PYTHIA 8 A2 tune 
describes the data at a similar level as the pythia 6 
AMBT2B tune. With the exception of tune Zl, all of 
the models significantly overestimate the fraction of dijet 
events. For the transverse sphericity, shown in Fig. 3(b), 
the different PYTHIA tune predictions, other than Zl, are 
similar. All models tend to reproduce the data selected 
with p!£ ad > 5 GeV better, as seen in Figs. 2(e), 2(f) and 
3(c). The Zl tune again provides the best description of 
the data, followed by Herwig++ and pythia 8 A2. For 
this distribution, the DW tune also provides the same 
level of agreement as the other tunes. 

The differences of the MC distributions with respect to 
the measured event shape distributions do not however 
yield more than a 5 — 10% shift in their corresponding 
mean values. With the exception of pythia 6 DW, the 
MC models seem to predict the plateau value reasonably 
well in Fig. 4. The mean values of all event shape variables 
increase up to A c h of about 30, or up to X^Pt of about 
30 GeV and tend to saturate at higher values of A c h and 

The ALICE collaboration has measured the transverse 
sphericity distribution in inelastic 7 TeV pp collisions [13], 
selecting charged particles with \r]\ < 0.8. The behavior 



of mean transverse sphericity as a function of multiplicity 
exhibits a similar behavior to that observed here, with 
the data lying at higher values than predicted by MC 
models. 



X. CONCLUSIONS 

The event shape observables, transverse thrust, trans- 
verse thrust minor, and transverse sphericity, have been 
measured requiring at least six charged particles per event, 
in low pile-up minimum bias proton-proton collisions at 
y/s = 7 TeV. The distributions and mean values have 
been compared to a number of MC model predictions 
for various tunes of the minimum-bias and underlying 
event models. The dependence of the event shapes on the 
number of charged particles, on the sum of pt and on 
the leading charged particle pr has been studied. 

The distribution of transverse sphericity shows a tran- 
sition from spherical events toward a dijet structure as 
p!p d increases, while transverse thrust and thrust minor 
are much less sensitive. The dependence of the event 
shape mean values as functions of and J^Pt is sim- 
ilar. For all variables, a slow variation is observed at 
low multiplicity and a saturation towards more spherical 
events at higher multiplicity. No evidence is observed for 
a dijet structure at the highest measured values of iV ch 
and J2pt- All MC generators underestimate the fraction 
of spherical events and none reproduces accurately the 
event shape distributions. The MC tunes based on the 
properties of the underlying event show in general better 
agreement with the data than those based on the inclu- 
sive distributions measured in minimum-bias events. The 
pythia 6 with the Zl tune provides the most accurate 
description for the distributions presented in this analysis, 
but the MC agreement could still be improved. These 
measurements provide information complementary to in- 
clusive particle distributions and thus they are useful for 
improving the MC description of minimum-bias collisions 
at the LHC. 



d 



6 Er- 

5 Er- 
4 — 



<s = 7 TeV ATLAS 

P T > 0.5 GeV, N > 6 and h| < 2.5 , p^"" > 0.5 GeV 

- Data 2010 

- PYTHIA 6 AMBT2B 
■■ PYTHIA 6 DW 

- PYTHIA 6 Z1 
■ PYTHIA S A2 

- Herwig++ UE7-2 




0.35 
T ch 



"O I -o 

I > 



a 

a 



6 — 

5 E- 

4 ^- 

3 ^~ 
2 ^- 
1 — 



is = 7 TeV ATLAS 

P T > 0.5 GeV, N h > 6 and h| < 2.5 , p*"" > 0.5 GeV 

— Data 2010 

-. PYTHIA 6 AMBT2B 

— PYTHIA 6 DW 

— PYTHIA 6 Z1 
- PYTHIA 8 A2 

— Herwig++ UE7-2 





(a) 



(b) 



T3 | <-> 

I > 

T— <0 



Q 

O 



7 — 

6 EH 




is = 7 TeV ATLAS 

p T > 0.5 GeV, N h > 6 and h| < 2.5 , p*"" > 2.5 GeV 

- Data 2010 

- PYTHIA 6 AMBT2B 
■■ PYTHIA 6 DW 

- PYTHIA 6 Z1 

- PYTHIA 8 A2 

- Herwig++ UE7-2 i" 



. 1 . 1 1 . . 


















0.05 


0.1 0.15 0.2 

(c) 


0.25 0.3 0.35 

xf 



mo* 

^ i- 

"O I "O 



6 — 

5 ^- 
4 ^- 
3 E~ 
2 r - 

1 — 



is = 7 TeV ATLAS 

p T > 0.5 GeV, N h > 6 and h| < 2.5 , p*" d > 2.5 GeV 
■ Data 2010 



PYTHIA 6 AMBT2B 
PYTHIA 6 DW 

- PYTHIA 6 Z1 

- PYTHIA 8 A2 
Herwig++ UE7-2 




(d) 



-a \T3 



a 

o 



7 Er 
6 E 
5 =- 
4 



Vs = 7TeV ATLAS 

p T > 0.5 GeV, N > 6 and h| < 2.5 , p*"" > 5 GeV 

-•- Data 2010 

PYTHIA 6 AMBT2B 
-— PYTHIA 6 DW 

PYTHIA 6 Z1 

PYTHIA 8 A2 

Herwig++ UE7-2 




nivT^ i i i i i i ■ i i i i i ■ ■ i i i ' ■ ■ ' ' ' ' ' ' ■ ' ' 


'_ t-'~- 






; 






0.05 0.1 0.15 0.2 


0.25 0.3 0.35 

xf 


(e) 





-o I -o 



Q 

O 




-pch 
' M 



(f) 



FIG. 2. Normalized distributions of transverse thrust (left) and transverse thrust minor (right) for particles with px > 0.5 GeV 
and \r)\ < 2.5 for different requirements on the transverse momentum of the leading charged particle, p^ d (top to bottom). The 
error bars show the statistical uncertainty while the shaded area shows the combined statistical and systematic uncertainty. 
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FIG. 3. Normalized distributions of transverse sphericity for 
particles with pt > 0.5 GeV and |tj| < 2.5 for different require- 
ments on the transverse momentum on the leading particle, 
p!}; ad (top to bottom). The error bars show the statistical un- 
certainty while the shaded area shows the combined statistical 
and systematic uncertainty. 
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FIG. 4. Mean values of transverse thrust, transverse thrust minor and transverse sphericity (top to bottom) with px > 0.5 GeV 
and 1 77 1 < 2.5 versus charged particle multiplicity of the event (left column) and versus charged particle transverse momentum 
scalar sum of the event (right column) for transverse momentum of the leading particle p^ ad °f 0-5 GeV. The error bars show 
the statistical uncertainty while the shaded area shows the combined statistical and systematic uncertainty. 
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